COMPUTERS 


Electro/Scientific Industries 


7524 S.W. MACADAM AVE.+ PORTLAND I9, OREGON JANUARY 1960 
REPLACES AUGUST 1959 


A NEW COMPUTER 
FOR 


ALGEBRAIC FUNCTIONS 
OF A 


COMPLEX VARIABLE 


MERLE L. MORGAN, Ph.D. 


Director of Research and Development 


Litho in U.S.A. 


7, _—— 
has =. eter — 


ravi a] a) 6 o> 


ABS TRACT 


A NEW COMPUTER 
FOR 
ALGEBRAIC FUNCTIONS OF A COMPLEX VARIABLE 


The ESIAC algebraic computer is presented as a new tool for 
use in control system analysis and design. The ESIAC is a special- 
ized analog computer which serves the designer as a '"'two-di- 
mensional slide rule" for functions of a complex variable and (heise 
pole-zero plots. 

The features of the ESIAC are described and a number of 
applications discussed, including its use in automatic plotting of 
root locus and frequency response plots, in the evaluation of res- 
idues, in factoring polynomials, in network synthesis, in the ana- 
lysis and design of carrier systems, systems with dead time and 
sampled data systems, in the treatment of non-linear systems by 
the describing-function method, and in the design of systems by 


power-density spectra methods. 
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A NEW COMPUTER FOR 
ALGEBRAIC FUNCTIONS OF A COMPLEX VARIABLE 
by Merle L. Morgan 


Introduction 


The problems associated with the development of modern 
guidance systems, process controls, and other servomecha- 
nisms have placed increasingly greater demands on the design 
engineer. The complex nature of these problems has resulted 
in the development of more effective methods of analysis and 
synthesis for their solution. For example, in the last decade 
there has been rapid development and extensive use of Laplace 
transform techniques. The representation of functions by pole- 
zero diagrams and root locus plots in the study of the transient 
and steady-state responses of systems has become widespread. 

The advantage of using Laplace transform techniques is that 
complicated integro-differential equations in the time domain are 
converted to simple algebraic equations in the frequency domain. 
This allows the engineer to do his mathematical manipulation in 
terms of algebraic functions rather than differential equations. 
These algebraic techniques are applicable not only to ordinary 
linear servomechanism design and network synthesis but also to 
the analysis and design of carrier systems, systems with dead 
time, sampled data systems, and various nonlinear systems by 
the describing function method, and to the design of systems by 


power density spectra methods. 
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Evaluation of Algebraic Functions 

The algebraic functions of a complex variable and associated 
pole-zero plots encountered in such problems are commonly eval- 
uated graphically. As the systems become more complicated, 
however, the graphical techniques become quite time consuming 
and the use of a computer is very desirable. 

Unfortunately, the design engineer has not previously had at 
his disposal a simple computer for algebraic functions described 
in terms of their zeros and poles. A large digital computer can 
be programmed to handle such functions but the use of its equip- 
ment and personnel is costly, and its operating time usually must 
be scheduled in advance. For these reasons , the designer cannot 
quickly obtain a solution for one trial set of parameters and then 
proceed immediately with a solution for a new set chosen upon 
inspection of the first solution. 

Analog methods are typically much simpler and more suit- 
able for calculations where the high accuracy of digital methods 
is not required. Conventional analog computers, however, are 
designed to solve differential equations involving functions of the 
real variable time and are not easily adapted to functions of a 
complex variable. 

Two-dimensional electric potential distribution methods on 
the other hand have a very simple, direct relationship to functions 


of a complex variable. Potential distribution analogs have been 
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used very effectively as an aid in visualizing functions in a 
complex plane. They offer the additional advantage that any con- 
formal transformation of the coordinate system can be made with- 
out affecting the function, providing proper boundary values are 
maintained. Their practical application to computation, however, 
has been limited in the past to the direct potential analog ofa 
function, in which current sources and sinks are placed at the 
zeros and poles of the function. This provides a direct analog of 
the logarithm of the magnitude of the function, but no similar 
analog can be constructed in general for the angle of the function. 

The new ESIAG’ Algebraic Computer” provides for the first 
time analog computation of both the magnitude and the angle of a 
function of a complex variable. This is accomplished by means 
of a combination of electric potential distribution and lumped 
circuit elements. With the ESIAC the advantages of analog com- 
putation are available for evaluation of the functions of a complex 
variable encountered in all the algebraic system design methods 
previously mentioned. In transient analysis problems, the ESIAC 
provides an improvement on the graphical root locus technique. 
The ESIAC also measures the residues at the first order poles 
of a function, thus providing the partial fraction expansion for 
the transient response of a system. For frequency response cal- 
culations, the frequency axis of the complex plane is scanned to 
plot the gain and phase shift of various portions of a system as a 


function of frequency. 


* Trade mark of Electro-Measurements, Inc., Portland, Oregon, 
We Ses 


° U.S. A. and foreign patents pending 
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Factored Algebraic Functions 


Algebraic functions of a complex variable are sometimes ex- 
pressed as the ratio of two polynomials. When the polynomials 
are factored, however, such a function can always be expressed 


in either of the forms: 


n n n Nn, 
be 0 | a 2 = 3 
n n 
= 0 i 
or 
n n n n 
0 s 1 s 2 Ss « 3 
F=+ K.s (1 - ——) (1 -——) (1 - ——) Seis 
2 8) 84 83 
n n 
0. s i 
st K,5 i (1 a (2) 
; 1 1 
in which 


F, s and the 8,'s represent complex vector quantities. 
K, or K,: oy and the n,'s represent scalar quantities 
which can be either positive or negative. 
Any rational algebraic function can be factored into either of 
the above forms, and the Laplace transform of any linear differ- 
ential equation can therefore be expressed in this form. The 


ESIAC handles functions expressed in the form of either of these 


two equations. 
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Such factored functions are evaluated by recognizing that the 
magnitude of F is equal to the product of the magnitudes of the 
factors on the right hand side of the equation (its log magnitude is 
equal to the sum of the log magnitudes on the right) and the angle 
of F is equal to the sum of the angles of the factors on the right 
hand side of the equation. 

In practical engineering problems, it is usually necessary to 
determine the values of IF 1 (magnitude of F) and ZF (angle of F) 
corresponding to each of a large number of values of s. Some- 
times the values are used directly as in plots of gain and phase 
as a function of frequency. In other cases, a large number of 
trial values of s must be substituted in the equation to find the 
values of s which satisfy the equation when the value of F is.spec- 
ified. These values are the roots of the equation. This technique 
is used in the analysis of control systems by the root locus method. 
The same approach is also applicable to the factoring of high order 
algebraic polynomials, an operation which is often required as a 
preliminary to the application of the root locus method to the 
equation of a system. 

The numerical and graphical methods for solving such prob- 
lems are tedious since the magnitude and angle of (s - s.) or 
(1 - —) must be determined for each factor for every trial value 
of s. The unique advantage of the ESIAC is that a change in s in- 


volves setting only the dials for the magnitude and angle of s. The 
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values of the s,'s and other parameters are established at the 
start of the problem and require no further attention as the values 
of IF | and ZF are determined at various values of s. 
Description of the ESIAC 

The ESIAC consists of a console with a plotting board at the 
center and panels containing the principal controls conveniently 
arranged around it , as shown in Figure 1. A close-up of the 
operating controls is shown in Figure 2, and in Figure 3 the 
covers at each side of the plotting board have been raised to 
show the two electric potential analog sheets and the frames 
for positioning and moving probes over them. The supply of 
probes and probe clamps in the drawer is also shown in Figure 3. 

The ESIAC contains two separate analog circuits--one for 
[#1 and the other for ZF. In setting up a problem, a pattern of 
probes is placed in the frame over each sheet at positions corres- 
ponding to the s,'s (the zeros and poles) in the equation. The 
frames are then moved over the sheets to positions corresponding 
to values of the complex variable s. The voltage distributions on 
the sheets are such that the voltage at each probe represents on 


one sheet the magnitude and on the other sheet the angle of one 


factor (s - 8.) or (1 - - ). The exponent n, on each factor is 
he i Py | 
established by selecting the value of an electric coupling element 


in the clamp which attaches the probe to the frame. Dials on the 


panel provide for setting the values of F, K ands, and for selecting 
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FIGURE 2. ESIAC OPERATING CONTROLS 


FIGURE 3. ANALOG SHEETS AND PROBES 
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the sign of K, the value of the exponent Ny and the form of the 
general factor (s - s;) or (1 - =). 

When the controls edoleas all the parameters in the 
equation are adjusted to any combination of values for which the 
two sides of the equation are equal, a pair of electrical null indi- 
cators, one for the magnitude and one for the angle, show the con- 
dition of equality. After the values of all given or assumed par- 
ameters have been set, the equation is solved simply by adjusting 
the value of the unknown until null balances are obtained. 

For root finding or root locus plotting, the motion of the probe 
assembly is directly coupled to the stylus on the plotting board. 
For automatic plotting, the probe assembly is moved to scan values 
of s, and whenever a magnitude or angle balance is passed, the 
stylus automatically marks the position on the plotting board. A 
line joining the marks for magnitude balance is the solution of the 
magnitude equation; the line for angle balance is the solution of 
the angle equation. The value of s at each intersection of these 
lines is therefore a root of the complex equation. 

To use the ESIAC for frequency response plotting, the po- 
sition of the probe assembly in the coordinate representing Zs 
is fixed at 90°, and its motion in the Is! direction is coupled to 
the plotting assembly in the horizontal direction. The vertical 
motion of the plotting stylus is coupled to the IF! dial for gain 


plotting or to the 7F dial for phase plotting. When the circuit is 
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balanced at various values of Isl, the stylus plots the log gain or 
phase angle versus log frequency. 

Since the plotting stylus is directly coupled to the dials, it is 
not necessary to read the dials while making these plots. The 
stylus is simply moved about over the plotting paper and the auto- 
matic marking circuit plots the curves. 

Logarithmic Coordinates 

In the ESIAC, each parameter is represented in terms of its 
magnitude and angle. The magnitude scales are logarithmic, sO 
that all values are represented with the same relative accuracy, 
as on a slide rule. The angle scales are linear. The complex 
plane on which values of s are plotted is therefore the log s 
plane rather than the s plane. 

There are both electrical and mathematical advantages in 


using the log s plane. Electrically, the use of the log s plane 


oo 


makes possible the inclusion of both the (s - s.) and (1 - 
forms for the factors in Equation 1 and Equation 2. The na 
equations differ only in the value of K, but it is convenient to be 
able to use either K which may be specified in a problem. 
Probably the most important advantage of the log s plane is 
that values of s over a wide range can be represented on the 


same plot. The ESIAC covers a range of four decades of s and 


eight decades of s.'s, thus including the effect of zeros and poles 
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two decades beyond the ends of the plot. The use of the log s plane 
results in uniform relative accuracy over the whole four decades, 
and in compatibility between the scales for root loci and the scales 
for Bode plots. 

In the log s plane, any complex quantity is repeated every mul- 
tiple of 360° in the Zs direction. Therefore only one strip of 360° 
width is required for plotting all complex quantities. Any multiple 
of 360° can be added to or subtracted from the angle used to specify 
a complex quantity without changing the quantity. When analog 
voltages appear in the process of computation which are outside 
the 360° range of the ZF scale, a stepping switch automatically 
adds or subtracts multiples of 360° until the voltage is brought 
within the range. 

Two types of graph paper, shown in Figure 4, are provided 
with the ESIAC for plotting values of s in the log s plane. The 
semi-logarithmic graph paper, ESIAC Form 46, has a logarithmic 
scale for Isl and a linear scale for Zs, corresponding directly 
to the scales on the ESIAC. These coordinates are the loga- 
rithmic transformation of the polar coordinates of the s plane. 
The transformation of the rectangular coordinates of the s plane, 
Re [s] and Im [s| , produces the curvilinear coordinates of 
ESIAC Form 47. The four arms of the axes in the s plane, +, 4), 


-~ and -j, become four parallel horizontal lines in the log s plane. 
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To read the value of s at any point on Form 47, it is necessary to 
follow the curves from the point to their intersections with the two 
adjacent parallel axes and read the coordinates there. Since all 
the axes have the same logarithmic scale, their values are marked 
only along the bottom and top edges of the graph paper. 

In pole-zero diagrams or root locus plots, poles lying on the 
180° axis represent critically damped responses. The damping of 
a response decreases as its poles move farther from 180° and 
reaches a state of continuous oscillation when the poles are at 
90° ana 270° respectively. Since the damping factor is the nega- 
tive of the cosine of Zs, contours of constant damping factor are 
parallel horizontal lines across the sheet. Both Forms 46 and 47 
are provided with damping factor scales at their right margins as 
a convenience. 

Frequency response plots on the ESIAC are simply plots of IF 1 
and ZF versus Isl (angular frequency w) or J (frequency f), when 
Zs =90°. These appear as the familiar Bode plots of logarithmic 
gain or linear phase versus a logarithmic scale of frequency or 
angular frequency. 

Typical ESIAC Applications 

The ESIAC is applicable to a wide variety of problems en- 

countered in control system analysis and design. It can be used 


on any problem which can be expressed in the form of equation 
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(1) or equation (2). Either F or s Say be the unknown in the 
equation. The use of the ESIAC in a number of typical eonteol 
system problems will be discussed briefly in the following 
sections. The discussions are not intended to be complete, but 
the general approach to various problems is indicated. Since 
the sane mathematical methods can all be found in standard 
control.system textbooks, no references will be given. 


Root Locus Plots 


. One of the most frequent applications of the ESIAC is in the 
construction of root locus plots. In general, a root locus plot is 
a euewe: of the locus of the roots of an algebraic equation as any 
parameter is varied. The term is most commonly applied to the 
determination of the poles of the closed-loop transfer function of 
a system by setting the denominator of the transfer function equal 
to zero and plotting the locus of the roots of this equation as the 
gain parameter K is varied. Since K affects only the magnitude 
balance of the equation, the locus is simply the solution of the 
angle equation, which is plotted directly on the ESIAC,. An ex- 
ample of such a root locus plot for a moderately complicated 
control system is shown in Figure 5. 

For a discussion of further details of such a problem, a much 

simpler example will be used. Its block diagrams and equations 
are given in Figure 6, Since two parameters K, and K. are to be 


chosen, two or more root locus plots will be made. The first, 
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3? 8, and 8, of 


the portion of the system including the forward path and the unity 


shown in Figure 7 , is used to locate the poles s 


feedback path. The variable parameter in this plot is K.- A 
second plot can then be made for any chosen value of K, to deter- 


mine the poles Sp: 8 and sg of the complete system as a function 


7 
of Kis For a choice of K. = 120, the second plot is shown in 
Figure 8 . Similar plots can be made if desired for other values 


s,ands,in Figure 7. 


of K, by starting from their values of S358, 5 


Note that the numerical values of s and Sp need not be read. 


ai aA 
Their graphical locations on Figure 7 are simply transferred 
directly to Figure 8 . Assuming that K, is to be chosen for the 
maximum damping factor of S 6 and 80, their locations are deter- 
mined by inspection in Figure 8. At the points indicated, the 
value K, = 0.18 is read on the K scale of the ESIAC. The 
position of Se for this value of K, is found to be as shown in the 


same figure. The numerical values of these roots can be read 


on Figure 8 as: 


So = =G.- J40 
Se = -5.4 


These are the final closed-loop poles for the system shown in 


Figure 6. 
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Block Diagram of System: 


Amplifier, Motor and Load 


Inpu K 


6. cs 


s 
eg! ame 


Position Feedback 


Rate Feedback 


Sy and S, are the roots of the equation 


in which $3, 
K 
1 


s s ne 
tania a Soe 


For example, when Ss) =- 6 ands, = - 12, the roots are plotted as a function of 


K, in Figure 7. 


Second Reduction of Block Diagram: 


1 


Output 
s s s 
Ea ay ae 


8 


in which So. 87 and Sg are the roots of the equation 
s 
K 


(1 - =-)(1- =-)(1 - =) 
“3 $4 55 


For example when K, = 120, the roots are plotted as a function of K, in Figure 8. 


For maximum damping, set K, = 0.18. This places s, and S7 at - 6 + j40 for a 


damping factor of 0.15, and places Sg at - 5.4. 


For these pole locations, the frequency response curves are plotted in Figures 9 


and 10. 


FIGURE 6. ANILLUSTRATIVE PROBLEM FOR THE ESIAC 
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Transient Response by Partial Fraction Expansion 


In making the inverse Laplace transformation of a 
response function by a partial fraction expansion, coefficients 
can be found by removing each pole factor and setting the 
variable to’ the value of the pole removed. The ESIAC 
performs these two operations simultaneously and gives the 
value of the coefficient or residue. The computer values can 
be used to find the complete partial fraction expansion even 
in the presence of high order poles. 


In the system of Figure 6 assume that the input is a 
unit step 05 == The output of the system is then 


9000 


oa (s - s,)(s - sy) (8 - 8) 8 
6 7 8 


or 


By setting the factored function 06 on the ESIAC, the residues at 


s. ands, can be read: 


moe 8 
A, = 0.07 280° at sp = - 6 + j40 
A, = 0.07 £-80° at 8s, =- 6 - j40 
Oo 
Ag = 12180 at Bg = 264 


By inspection of the function: 


A rl at s =0 
e) 
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Taking the inverse Laplace transform term by term, the output as 
a function of time is 


oO 
0 (t) = eee fe ,(=)] = 0.1467 aon 


6t 
27 ) 


cos (80° + 40t eg ey 


Frequency Response Plots 


If the frequency response of the system to a sinusoidal input 
is desired, it is easily plotted on the ESIAC. The final closed- 
loop gain and phase of the system of Figure 6 are simply the mag- 


nitude and angle respectively of ‘the transfer function 


ea Fa oe ey | caer 8 
i m = 5 


in which s = jw, and S, 84 and 8, are as shown in Figure 8. 


In terms of the ESIAC controls, s = jw means that Zs = 90° 


Is] 
ot 


are indicated directly on the ESIAC scale, and graph 


= f (frequency). Both 


and that Isl = w (angular frequency) or 


Is] 
aT 


papers are provided for use with the ESIAC in plotting IF 1 andZF 


Ist and 


as a function of either w or f. The plotting stylus in the ESIAC 


i 


can be directly coupled to IF | or ZF for plotting frequency response 


plots. The resulting curves as a function of f for this example 


are shown in Figures 9 and 10. 


Factoring Polyn omials 


Problems are frequently encountered in which polynomials 
must be factored. Any polynomial can be expressed in either of 


the factored forms: 
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the s;'S are S&, 
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, s. 

i =f 1 


in which the s.'s are the roots of the equation formed by setting 
the polynomial equal to zero. 

This factoring can be accomplished on the ESIAC bya 
sequence of root locus plots, in which each plot is used to re- 
place a group of terms in the polynomial by its factored equiv- 
alent. For example: 


2 
a, t 4,8 +ai,sS ta sees se pe s° + a,s cate 


l 2 3 4 5 
. Ss Ss Ss 4 5 6 
See eee any ign eg + as + ays ee 
6 
za. (1 -—=-)(1 - =—)(1 - S—)(1 - =) - 4+ a,s°+... 
Oo 84 S. 86 s_ 89 6 


Sal 
—- eo 


in which s,,; 8 are the roots of the equation 


l 3 ands 


S| 


sta go 45 s? = 0 


lly 2 3 


S4? Sp, 86s 82 and S, are the roots of the equation 


Ss Ss Ss 4 5 
ao (1 ore 0. - le popeer a,S +a,s° = 0 
et cetera. 
The equation for s, through s, is solved on the ESIAC by a root 


1 3 


locus plot after rewriting it 
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The equation for s, through s, is similarly solved after rewriting 


4. 8 


it 


a, (1 a - = - =) 


3 
a_s 
as” 1+ 
4 


This procedure is repeated until the whole polynomial is factored. 


Network Calculations 

In network calculations, a variety of problems arise which 
require factoring or other operations involving root finding by 
locus tracing techniques. For example, when two factored 
transfer functions are to be added and the result expressed in 
factored form, the numerator factors may be obtained by a root 
locus plot. In this nvebicm: 


A(s) 
Bis) 


i" C(s) _ A(s) D(s) + B(s) C(s) 
Dis) - B(s) Dis 


The denominator is already factored but the factors of the numer- 
ator must be found by solving the equation 

A(s) D(s) + B(s) C(s) = 0 
This equation can then be expressed in the form 


A(s) D(s) _ 


‘B(s) C(s) = 


and the roots found on the ESIAC. 
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Another similar problem is the determination of loading 
effect. When a network function Z(s) for infinite impedance out- 
put is known and a specified load resistance R(s) is to be connected, 
the transfer function of the cascaded system requires solution of 
the equation 

1 + ee = 0 
in which Z(s) is the driving point impedance of the unterminated 
network and R(s) is the terminating impedance. 
Representing Dead Time on the ESIAC 

In a system with dead time, such as a transport lag or a delay 
line, the expression erie appears. For calculations on the ESIAC, 
this expression is approximated by II (1 - Seay in which the 
zeros and poles Ss, are located ata Peon above the 
maximum signal frequency. This approach is the same of that of 
approximating a distributed delay line or filter by a network of 
lumped circuit elements. 

Several approximations of this form may be made with 
different locations and orders of zeros and poles. An example 
using four zeros and four poles is given in Figure ll. 

Desc ribing Functions for Non-Linear Systems 

Many systems containing a non-linear element can be analyzed 

by the describing function method. In this method a Fourier 


analysis for the output of the non-linear element with a sinus- 


oidal input is determined as a function of signal amplitude. It 
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18 
is assumed that only the fundamental component of the output of 
the non-linear element contributes to its input. This assumption 
is reasonable in many systems as a result of harmonic attenuation 
in the rest of the loop. 

The gain and phase shift of the non-linear element at the fun- 
damental frequency can be combined with the rest of the system in 
the same manner as for a linear element except that the transfer 
function for the non-linear element varies with amplitude. Asa 
result, the analysis must be repeated for signals of various am- 
plitude. The stability of a feedback loop containing such a non- 
linear element is determined by calculating first the sinusoidal 
transfer function for the linear part of the loop and then combining 
this result with the gain and phase shift of the non-linear element 
at various amplitudes. This analysis can be made using any of 
the usual techniques--Nyquist diagram, Nichols chart, Bode plot 
or root locus plot. 

If the describing function contains a phase angle, the root 
locus method requires the plotting of loci for the linear portion 
of the transfer function not only for a 180° angle but also for other 
angles differing from 180° by the phase shift in the non-linear ele- 
ment. This added difficulty of construction has discouraged the 
use of root locus methods in such problems. With the ESIAC, 
however’, the method is much more attractive since it is a very 


simple matter on the ESIAC to plot a family of loci. 
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Carrier Systems 

In a carrier system, the closed-loop poles and system sta- 
bility can be determined from the transfer functions before and 
after sinusoidal multipliers for modulators and demodulators. 
The transform of the modulator output is the input transform 
shifted in frequency by the carrier frequency. This is a sum of 
functions whose zeros can be determined on the ESIAC. When 
the poles and zeros of the carrier frequency components are 
superimposed, the operation is repeated for the demodulator. 


Sampled Data Systems 

The analysis of sampled data systems requires a transfor- 
mation to the variable z = arP . For calculations on the ESIAC, 
its variable s will be considered as z rather than s. In the z 
plane the frequency axis is the unit circle Iz] = 1.0, not jw as 
in the s plane. The frequency response is the magnitude and 
angle along this line of constant Izl. In the logarithmic repre- 
sentation of the ESIAC, this contour is a straight vertical line 
at Izlel. 


Sampled data system functions are sometimes converted to 


functions of a new variable w defined by the relation 


in which z = are as in the previous method. These functions of 
w can be computed and plotted on the ESIAC in the same manner 


as functions of s. 
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Power-Density Spectra Problems 


In power-density spectra methods, the poles and zeros ofa 
power-density spectrum in the complex plane are given. There 
will be zeros and poles in all four quadrants of the complex plane. 
The spectrum can be measured along the frequency axis in the 
same manner as ordinary frequency response curves. 

If a given input spectrum is passed through a system or a 
filter, the output spectrum is obtained by superimposing the pole- 
zero pattern of the transfer function and its mirror image in the 
frequency axis on the pole-zero pattern of the spectrum. For the 
cross-power spectrum between the system output and another 
signal, the pole-zero pattern of the complex conjugate of the 
transfer function is added to the pole-zero pattern of the cross- 
power spectrum between the system input and the other signal. 
The output cross-power spectrum will have the phase as meas- 
ured along the frequency axis. 

When the spectra of signal and equivalent noise at a system 
input are given, together with the transfer functions of an unalter- 
able system components, the open-loop transfer function of the 
optimum control system can be determined by the Wiener-Hopf 
method. 


It is possible to determine the open-loop and closed-loop 


transfer functions of an operating process by measuring the spectra 
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21 
at the input and the output of the controller for two different con- 
troller gain settings. The transfer functions can then be obtained 
by making root locus plots from this data. 

All of the above power-density spectra methods involve 
algebraic operations which are conveniently performed on the 
ESIAC. 

Conclusion 

The ESIAC algebraic computer has been presented as a new 
tool for use in control system analysis and design. The ESIAC 
serves the designer as a "two-dimensional slide rule'' for functions 
of a complex variable and their pole-zero plots. The ease with 
which it locates roots and plots graphs makes possible rapid com- 
parison of various proposed stabilizing techniques and other design 
changes in terms of their algebraic transfer functions.. The ESIAC 
is particularly useful in a sequence of calculations where each sol- 
ution gives the designer further insight into the nature of the prob- 
lem and thus influences the direction which the investigation should 
take. 

It is hoped that this presentation of the features of the ESIAC 
and some of its applications will encourage greater use of these 
techniques, and perhaps even stimulate the development of better 


techniques than are now known. 
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